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We show that the motion of a laser-driven Bose-Einstein condensate in a high-finesse optical 
cavity realizes the spin-boson Dicke-model. The quantum phase transition of the Dicke-model 
from the normal to the superradiant phase corresponds to the self-organization of atoms from 
the homogeneous into a periodically patterned distribution above a critical driving strength. The 
fragility of the ground state due to photon measurement induced back action is calculated. 
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A thermal cloud of cold atoms interacting with a sin- 
gle mode of a high-finesse optical cavity can undergo a 
phase transition when tuning the power of a laser field 
which illuminates the atoms from a direction perpendic- 
ular to the cavity axis [TJ [2j [3l H]. Below a threshold 
power, the thermal fluctuations stabilize the homoge- 
neous distribution of the cloud, and photons scattered 
by the atoms into the cavity destructively interfere, ren- 
dering the mean optical field to be zero. Above threshold, 
the atoms self-organize into a wavelength-periodic crys- 
talline order bound by the radiation field which, in this 
case, is composed of the constructive interference of pho- 
tons scattered off the atoms from the laser into the cavity. 
The same phase transition can happen for Bose-Einstein 
condensed ultra-cold atoms, that is exempt from thermal 
fluctuations. For low pump power at zero temperature, 
the homogeneous phase is stabilized by the kinetic energy 
and the atom-atom collisions, a sharp transition thresh- 
old is thus expected [U [3]. In both examples the self- 
organization is a non-equilibrium phase transition with 
the distinct phases being stationary states of the driven- 
damped dynamics. 

In this paper we show that the Hamiltonian under- 
lying the spatial self-organization is analogous to the 
Dicke-type Hamiltonian [7\ and the transition to the self- 
organized phase can thus be identified with the superra- 
diant quantum phase transition [5]. Hence, the quantum 
motion of ultracold atoms in a cavity effectively realizes 
the Dicke model and may lead to the first experimental 
studies on this paradigmatic system. The accessibility of 
such a Hamiltonian dynamics is limited by the coupling 
to the environment. We explore how quantum noise infil- 
trates and depletes the ground state [H), imposing thereby 
a condition on the time duration allowed for the adiabatic 
variation of the macroscopically populated ground state 
by means of tuning an external parameter. 

We consider a zero-temperature Bose-Einstein conden- 
sate of a number of N atoms of mass m which is in- 
side a high- Q optical cavity with a single quasi-resonant 
mode of frequency u>c- Such a system has been re- 
alized and manipulated in several recent experiments 



PH HU [H EH HI E3- The atoms are coherently driven 
from the side by a pump laser field. The pump laser 
frequency u> is detuned far below the atomic resonance 
frequency uja, so that the atom-pump (red) detuning 
Aa — <jl> — loa far exceeds the rate of spontaneous emis- 
sion. One can then adiabatically eliminate the excited 
atomic level and the atom acts merely as a phase-shifter 
on the field. The dispersive atom- field interaction has 
a strength Uq = gfi/AA, where go is the single- photon 
Rabi frequency at the antinode of the cavity mode. We 
describe the condensate dynamics in one dimension along 
the cavity axis x, where the cavity mode function is 
cos kx. The motion perpendicular to the cavity axis re- 
quires a trivial generalization of the theory, and with a 
standing-wave side pump the self-organization effect oc- 
curs quite similarly in two-, and three dimensions [2 . 

The many-particle Hamilton operator in a frame ro- 
tating at the pump frequency ui and with H = 1 reads 
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where ^(x) and a are the annihilation operators of the 
atom field and the cavity mode, respectively. The cavity 
length is L, the detuning A c — w — u>c is the effective 
photon energy in the cavity. Atom-atom s-wave scatter- 
ing is neglected. Besides the dispersive interaction term 
Uq cos 2 kx, there is another sinusoidal atom-photon cou- 
pling term describing an effective cavity-pump with the 
amplitude rj t — flg /AA, where f2 is the Rabi frequency 
of the coupling to the transverse driving field. 

Self-organization is a transition from the homogeneous 
to a A-periodic distribution. The minimum Hilbert-space 
for the atom field required to describe this transition is 
spanned by two Fourier- modes, 
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where Co and c\ are bosonic annihilation operators. In 
the low excitation regime these two modes can be as- 
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sumed to form a closed subspace, so cjc + cjci = TV 
is a constant of motion giving the number of parti- 
cles. On invoking the Schwinger-representation in terms 
of the spin S with components S x — |(c{co + CqCi), 
S y = ^(c\cq — CqCi) and the population difference S z = 
|(c{ci — CqCo), the Hamiltonian Eq. (JTJ) confined into the 
two-mode subspace reads 

H = —5c a 1 a + lorS z + iy(a* — a)S x /VN 

+ ua^a^ + S z /N^j , (3) 

where Sc — — 2u, lur = hk 2 /2m, u = NUq/A, and 
y = \/2N f] t . In the first line one can recognize the fa- 
mous Dicke-model Hamiltonian with a coupling constant 
y tunable via the transverse driving amplitude r\ t - The 
last term is inherent to the BEC-cavity system, how- 
ever, it does not essentially change the conclusions to 
be drawn here as long as |m| < |<5 C |. This condition has 
to be anyway fulfilled so that the neglect of the next ex- 
cited Fourier-mode C2COs2/cx be justified in Eq. In 
the following, we will restrict the discussion to the pa- 
rameter regime which is needed for the self-organization 
[5]. That is, Sc < is required in order to avoid a dy- 
namical instability of the system, which arises from the 
motional heating induced by the slightly retarded cavity 
field [16]. 

The thermodynamic limit is defined as N — > oo and 
V — > oo, while the atom density p oc N/V is kept con- 
stant. The coupling constants u and y have been in- 
troduced such that they are proportional to the atom 
density, u cx N/V and y cx y/ N/V (there is a filling 
factor coefficient), and thus remain constant in the ther- 
modynamic limit. The ground state can be determined 
as in Ref. 8 . Let's use the Holstein-Primakoff repre- 
sentation in which the spin-N/2 degree of freedom is ex- 
pressed in terms of the bosonic operator b such that 5_ = 
VN-tfbb, S + = &t y/N - W6, and S z = Wb-N/2. The 
Hamiltonian transforms into 



H = -Sc a) a + ojRb^b + ua)ab^b/N 
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Next, let's employ the similarity transformation 
i)- 1 (P)D~ 1 (a)HD(a)D(l3), with the displacement op- 
erators, D(a) = explaa 1, — a*a} and D(P) = exp{/3tf — 
f3*b}, which does not change the spectrum of the Hamil- 
tonian. Formally, the transformation amounts to replac- 
ing b —>■ b + (3, a — ► a + a, and analogously for the hermi- 
tian adjoint operators in Q. The resulting Hamiltonian 
is then expanded up to second-order in the boson oper- 
ators. Note that the expansion of the nonlinear square 
root term can be performed only approximately, because 
the physically sensible Hilbert space for the operator b is 



truncated as b^b < N. Therefore the forthcoming results 
are exact up to 1/N. 

There is a pair of real Pa and «o such that the linear 
terms in the Hamiltonian in the displaced phase space 
vanish for a = iy/~N aa and (3 — \/N 0q. They obey 
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The trivial solution ao — Pa — always satisfies these 
equations, which corresponds to the physical state of a 
homogeneous condensate and no photon in the cavity. 
When a ^ and Pa ^ 0, the product of the two equa- 
tions leads to a second-order algebraic equation, 



Sc^r + y 1 
ulor + y 2 



0, 



(6) 



where we used that Sc (ulor + y 2 ) ^ 0. There is a phys- 
ically sensible solution in the range < /3q < 1 if and 
only if y > y crit = \Z-S c lor. Then, 
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For u = 0, which amounts to the normal Dicke model, 

2_ 2 

the solution is P 2 = - 2 y" A w ith the same critical value 
2/crit of pump amplitude. The light shift term does not in- 
fluence the threshold, because the zero mean fields make 
this term vanish below threshold. Note also that this re- 
sult for j/ cr ; t corresponds to the one calculated from the 
instability of the Gross-Pitaevski equation (GPE) [5], if 
this latter is taken in the g c — > (collisionless atoms) 
and k — > (no cavity loss) limit. The approach based 
on the GPE was exempt from the two-mode approxima- 
tion. On the other hand, owing to the simplicity of the 
two-mode model, the exact ground state, which might 
include macroscopic quantum correlations, can be well 
approximated from the Hamiltonian obtained up to the 
quadratic order: 

H = E + M a ] a + Jf '+ M » fet b 

+ (V 2 + ^ + iM c (ot-o)(6t + 6) , (8) 



where M Q = — Sc + up 2 
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M y = lor + ual- yaaPa — 1/2 , (9c) 
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FIG. 1: (Color online) Photon (dashed red lines) and motion- 
ally excited atom (solid blue lines) numbers in the ground 
state. Thick lines represent the contributions form the mean 
fields («q and 0o, these are the photon and atom excitation 
numbers divided by N, respectively), which can be the or- 
der parameters of the phase transition. Thin lines represent 
the incoherent excitations due to the squeezing, given by the 
quantum averages (a^a) and (b^b) taken in the ground state. 
Parameters: So = — WOur, u = — Q.Iur. 



The ground state is the vacuum state of the normal mode 
oscillators which have the cigcnfrcqucncies 

^2 = Mi+M^ ± ^ (Mg-M.M,)' ~ M()MyM 2 , (1Q) 

Below threshold, the energy gap to the first excited state 
vanishes as y/l — y] y cvl t on approaching the critical point 
(the exponent is thus 1/2). The ground state contains ex- 
cited Fock states of the uncoupled photon a and atomic 
b modes because of the squeezing (with the coefficient 
M x — My) and two- mode squeezing (with the coefficient 
M c ) terms. Below threshold M x — M y = 0, thus the 
ground state is simply the two-mode squeezed vacuum of 
the a and b modes, which is an entangled state [T7l 1 1 8] . 
The quantum average of the photon and atom motion ex- 
citations in the ground state can be seen in Fig.[T]togethcr 
with the mean field populations. The squeezing leads to 
a singular state at the critical point which amounts to a 
divergence of the incoherent excitations (a>a) and (tfb). 
Therefore, the detection of the continuous transition of 
the mean field amplitudes ao and /3q requires, e.g., ho- 
modyne light measurement or the observation of the in- 
terference between the homogeneous and the sinusoidal 
components of the atom field. 

The quantum phase transition associated with the 
ground state of the Dicke-type Hamiltonian must be in- 
fluenced by the cavity loss. The coupling to the en- 
vironment amounts to a quantum measurement of the 
coupled BEC-cavity system [9l [19], and has a back ac- 
tion on its state. Therefore, even at zero temperature, 
the ground state is being depleted, which process can 
be modeled as a diffusion. We calculate the rate of dif- 
fusion out of the ground state in the following. For a 



compact notation the variables are arranged in a vector 
R = [a, a) ib, w\. The Heisenberg equations of motion 
originating from the quadratic Hamiltonian ^ are lin- 
ear and are driven by quantum noise terms associated 
with the photon field decay, ^R = M.R + | ,where the 
matrix M contains the coupling between the bosonic cre- 
ation and annihilation operators, and the noise source is 
£ = [£, gi, 0, 0]. The only non- vanishing noise correlation 
function is = 2K<J(i-i'), where 2 K is the pho- 

ton loss rate. We neglect the dissipative — kci and — na) 
terms, because we are interested in the transient dynam- 
ics and not in the stationary regime of the system. Ini- 
tially, the dominant effect in irreversibly escaping from 
the ground state can be attributed to the infiltration of 
quantum noise (a diffusion process). 

The left and right eigenvectors lS k ' and A k \ respec- 
tively, of M can be used to expand the fluctuation vector 
R in terms of normal modes: R = ^ZuPkT^- By use 
of the orthogonality of the left and right eigenvectors, 
(Z , r^) = 5k. i, where (a, 6) is the scalar product, the 
normal mode amplitudes are obtained as pk = (^ fc \J?). 
They evolve independently as 



Pk(t) = e-^MO) 



/ e-^ t - t '^Q k (t')dt' , (11) 
Jo 

where the projected noise is Qk = (l^ k \i)- In the present 
Hamiltonian problem, the normal modes form hermi- 
tian adjoint pairs p + ,p i + with eigenfrequencies ±cj + , and 
P-,p^_ with frequencies from Eq. (10), respectively, 
where each pair corresponds to one of the normal mode 
oscillator of the quadratic Hamiltonian in ([8|. Second 
order correlations evolve as 



2k . —l[ k) if 



i(uj k +LJ[) 



(12) 



The first term represents the initial condition. The dif- 
fusion is due to the second term in which the linear time 
dependence can be written as being proportional to the 
sin(a;)/a; function, where x — {uik + u)i)8t/2. 

The total population in the excited states ab ove the 



ground state is given by (p + p+ + p_p-). Using (12 1 for 



its time evolution, the first term vanishes in the ground 
state, and, in the second term, LUk + loi = for both 
p\_p+ and p_p- terms. Thus the time evolution leads 
exactly to a linear increase of the excited population, 
the corresponding diffusion rate is plotted in Fig. [2] with 
dashed line. The singularity at the critical point reflects 
that the excitation energy of one of the normal modes 
tends to zero. 

Let us calculate the diffusion in terms of measurable 
quantities, such as the number of incoherent photons and 
motionally excited atoms, SN — (aJa + tfb). The inco- 
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FIG. 2: Diffusion out from the ground state. The rate of 
increase of normal mode excitations (dashed line) and that of 
photons and motionally excited atoms (solid line) with coarse 
graining |fc| _1 <C St <C Wr 1 . The almost overlapping dashed- 
dotted line is derived from an adiabatic elimination method 
and is given by the analytic result of Eq. (151. Parameters: 
8c = — lOOwfl, u = —O.Iljr. 



herent population evolves as 

SN(t) = J2{Pk(t)pi(t)) (ri k \{ 



+ r (*L(i) 
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By using Eq. (12) and by approximating the sin(a;)/a: 
function, the diffusion rate becomes 
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where is the Heavyside-function. If the "time step" 
St is shorter than any of the time periods k>± , none of 
the pairs (k,l) is cut off by the Heavyside-function in the 
sum ( 14 1 . Then, it follows from the completeness relation 

that the depletion rate is zero. On 



k' i L j ~ °y 
such a short time the quantum noise is associated with 

the photon field amplitudes a and a\ and normal-order 

products vanish at zero temperature. Diffusion in the 

populations is obtained when a coarse graining of the 

dynamics over a longer St is performed. We consider 

only the special case \Sc\ 3> ojr, when there is a large 

difference between the eigenfrequencies to±. The time 

step can be set such that |#c| _1 <C St <C , and the two 

pairs (k, I) with tok = = ±w_ ~ ±wr also contribute 

to the double sum in Eq. (14), in addition to the (k,l) 

pairs with uji — —tOk- Then, the departure from the 

ground state appears as a regular diffusion process in the 

motional excitation Fock space with a finite rate even at 

the critical point, which is plotted in Fig. [2] with solid 

line. 

As Sc is the far highest frequency in this example, 
there is an alternative avenue to the depletion rate which 



relics on the adiabatic elimination of the photon field 
variables, a and a\ following the method of Ref. [19]. It 
leads to an analytical approximation, which is plotted in 
dashed-dotted line in Fig. [2] 
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St 
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Below threshold, the diffusion rate is about 
(k/\Sc\) (y/z/crit) 2 - Adiabatic following of the 
ground state by means of slow variation of the detuning 
(or the pump amplitude) requires that the smaller 
excitation frequency, w_ ujr^/1 — (y/j/ C rit) 2 > be much 
larger than the diffusion. The use of large detuning 
\Sc\ 3> K removes the time limitation imposed by the 
quantum noise. Although the critical point can be 
adiabatically approached only as close as y/y C rit -C 1, 
which is a generic feature of criticality, it is an intriguing 
possibility that matter wave and light field entanglement 
can be adiabatically created by making use of this 
critical system. 

In conclusion, we have shown that the zero temper- 
ature limit of the atomic self-organization in a cav- 
ity corresponds to the quantum phase transition given 
by the Dicke model. This connection is principally 
different from the proposals where some internal elec- 
tronic dynamics of the atoms in the cloud is involved 
HQ HH [23] ■ The key point here is that the ener- 
gies of the decoupled systems are much lower, the atom 
field excitation is being in the recoil frequency range of 
kHz, that of the photon field is broadly tunable, and the 
critical regime can be addressed in currently running ex- 
periments [2"4"] . 
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